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ABSTRACT 

We combine SDSS and WISE photometry for the full SDSS spectroscopic galaxy sample, creat¬ 
ing SEDs that cover A=0.4-22/im for an unprecedented large and comprehensive sample of 858,365 
present-epoch galaxies. Using MAGPHYS we then model simultaneously and consistently both the 
attenuated stellar SED and the dust emission at 12/rm and 22^tm, producing robust new calibrations 
for monochromatic mid-IR star formation rate proxies. These modeling results provide the first mid- 
IR-based view of the bi-modality in star formation activity among galaxies, exhibiting the sequence 
of star-forming galaxies (“main sequence”) with a slope of dlogSFR/dlogM * = 0.80 and a scatter 
of 0.39 dex. We find that these new star-formation rates along the SF main sequence are systemat¬ 
ically lower by a factor of 1.4 than those derived from optical spectroscopy. We show that for most 
present-day galaxies the 0.4-22/im SED fits can exquisitely predict the fluxes measured by Herschel at 
much longer wavelengths. Our analysis also illustrates that the majorities of stars in the present-day 
universe is formed in luminous galaxies (~ L*) in and around the ‘green valley’ of the color-luminosity 
plane. We make the matched photometry catalog and SED modeling results publicly available. 
Subject headings: catalogs — galaxies: statistics — galaxies: stellar content — galaxies: star formation 
— infrared: galaxies — 


1. INTRODUCTION 

Some of the most basic and important constraints on 
galaxy formation models are the present-day stellar mass 
function and the distribution of star formation among 
galaxies with different masses. The Sloan Digitital Sky 
Survey (SDSS, Blanton et al. 2005; Adelman-McCarthy 
et al. 2008; Padmanabhan et al. 2008) has provided the 
measurements that underlie our current knowledge of the 
stellar masses (M*) and star formation rates (SFRs) of 
large samples of galaxies (e.g., Brinchmann et al. 2004). 
The SDSS multi-wavelength ( ugriz ) imaging has been 
used to estimate luminosities and mass-to-light ratios 
(e.g., Baldry et al. 2012). The resulting stellar masses 
have been demonstrated to correlate tightly with total, 
dynamical mass estimates (Taylor et al. 2011) with a 
scatter of only 0.13 dex. Brinchmann et al. (2004) use 
the Chariot & Longhetti (2001) photoionization model 
to convert nebular emission line fluxes from SDSS spec¬ 
troscopy into SFRs. Salim et al. (2007) provided SFRs 
based on GALEX UV fluxes. But these previous SDSS 
studies did not account for the extra information enclosed 
in dust emission, when estimating both star formation 
rate and dust attenuation. Here we present alternative 
M* and SFR estimates by extending the photometric 
wavelength coverage of the SDSS spectroscopic sample to 
the near- and mid-infrared as enabled by WISE (Wright 
et al. 2010), which provides allsky coverage at 3 — 22/nn. 

The 3.4/mi and 4.5/rm bands sample the Rayleigh- 
Jeans tail of the stellar spectral energy distribution 
(SED), which avoids the contribution from hot, young 
stars which can dominate at shorter wavelenghts. In ad¬ 
dition, extinction is usually negligible in these bands. As 
a result, near-infrared luminosities provides fairly accu¬ 
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rate and precise stellar mass estimates, even in the ab¬ 
sence of any other photometric or spectroscopic infor¬ 
mation (Meidt et al. 2014). Many authors have used 
3 — 5/zm-photometry from Spitzer or WISE to compare 
with stellar mass estimates derived from photometry at 
2/im and below (e.g., Li et al. 2007; Zhu et al. 2010; Jar- 
rett et al. 2013). Here we combine SDSS and WISE data 
to generalize the use 3 — 5/mi photometric information 
in estimating stellar masses. 

Mid-infrared emission traces star-formation activity 
through the well-known correlation with PAH emission, 
sampled by the WISE 12/tm band, and through the 
correlation with thermal radiation from dust, sampled 
by the WISE 22/im band. A number of authors have 
compared star-formation rates from Brinchmann et al. 
(2004), based on optical emission lines, with mid-IR lu¬ 
minosities from WISE (e.g., Donoso et al. 2012; Lee et al. 
2013; Wen et al. 2013, 2014); these studies do not model 
the mid-IR luminosity to obtain SFRs that are indepen¬ 
dent of the emission-line based SFRs. 

Here we take the important step to include the full 
WISE photometry and employ SED modeling that con¬ 
sistently treats the stellar emission along with the dust 
extinction and emission. This will result in more ro¬ 
bust masses and star-formation rates for dusty galaxies 
and, in general, an alternative to the emission-line based 
SFRs from Brinchmann et al. (2004) for the full SDSS 
spectroscopic galaxy sample. This is not been attempted 
before for large samples of galaxies. Jarrett et al. (2013) 
provided a detailed multi-wavelength study of a small 
number of objects and will provide an important bench¬ 
mark to test our results. Brown et al. (2014) performed 
SED fits across the UV to mid-IR wavelength in an exer¬ 
cise similar in approach as what we present here, but for 
those authors the focus lay on producing a set of repre¬ 
sentative templates: stellar mass and SFR estimates are 
not presented or discussed. Cluver et al. (2014) analyze 
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Fig. 1.— Redshift distribution of the SDSS galaxy sample. 
98.24% of the objects are detected in at least one WISE band. The 
subsamples with > 2a detections in W3 (12/mi) and W4 (22/im) 
are highlighted. 

the mid-IR properties of the large GAMA sample, but 
there Ha-based SFRs were used as a calibrator for the 
mid-IR luminosities from WISE. As such, our study is 
the first to use the mid-IR luminosities of a large sample 
of galaxies (the full SDSS spectroscopic galaxy sample) 
to estimate their SFRs in a manner that is entirely in¬ 
dependent of optical emission line luminosities and other 
external calibrators. 

This paper is organized as follows. In Section 2 we 
will describe the SDSS and WISE photometric datasets, 
with a particular focus on total flux measurements for 
extended sources. In Section 3 we will describe MAG- 
PHYS (da Cunha et al. 2008), the SED fitting code, and 
test the robustness of the fitting results. In particular, 
we will test how well MAGPHYS can predict, with the 
available wavelength coverage of 0.4 — 22/im, the total 
IR luminosity. Furthermore, we will compare our M* 
and SFR estimates with the Brinchmann et al. (2004) 
measurements. In Section 4 we will describe the pub¬ 
licly released catalog with the fitting results for the en¬ 
tire SDSS spectroscopic sample (~800,000 galaxies). In 
Section 5 we will show, as an illustration, how star for¬ 
mation is distributed over galaxies with different masses 
and colors. 

We use AB magnitudes and adopt the cosmological 
parameters (Dm,Da,Ii)=(0.30,0.70,0.70) and adopt the 
Chabrier (2003) stellar initial mass function. 

2. DATA 

We use the SDSS spectroscopic galaxy sample as com¬ 
piled in the New York University Value-Added Galaxy 
Catalog (NYU-VAGC, Blanton et al. 2005; Adelman- 
McCarthy et al. 2008; Padmanabhan et al. 2008) 3 , 
which contains 858,365 galaxies with reliable redshift 
measurements (0BJTYPE=GALAXY, Zwarning=0), dis¬ 
tributed as shown in Figure 1. We adopt M0DELFLUX 
and M0DELFLUX_IVAR as the flux measurements and their 
uncertainties, with galactic extinction corrections as pre¬ 
scribed by Schlafly & Finkbeiner (2011). To account for 
the uncertainty in the extinction law, we propagate an 
uncertainty in Ry = 3.1 ± 0.2 in our flux uncertainty 
estimates. Minimum uncertainties are added in quadra¬ 
ture: 0.05 mag for the u band and 0.02 mag for the griz 
bands. This is aimed at preventing from small systematic 

3 http://sdss.physics.nyu.edu/vagc/ 


uncertainties to dominate the our goodness-of-fit assess¬ 
ment and the uncertainties in our derived parameters. In 
case of catastrophic failures in the photometric measure¬ 
ments, which occasionally occurs in the u and z bands, 
we omit these data points from the fits. 

The SDSS parent sample is cross-matched with the 
A11WISE source catalog 4 . The vast majority of galax¬ 
ies have counterparts (98.24%), identified as the bright¬ 
est object within a search radius of 6 arcsec, similar 
to the WISE PSF. In line with the approach taken by 
Donoso et al. (2012) and Yan et al. (2013), most matches 
(99.08%) are found within a radius of 3 arcsec and the 
majority (91.20%) do not have multiple matches within 
6 arcsec. We provide a flag (FLAGJV) that indicates 
if there are one or more matches, but we do not exclude 
any matched sources a priori. Various flux measurements 
are provided by the A11WISE catalog. In order to min¬ 
imize the effects of source blending we use the W?mpro 
and W?sigmpro flux measurements and their uncertain¬ 
ties. These are profile-fitted photometry measurements, 
performed simultaneously on neighboring sources, using 
the point spread function (PSF) as the source model. 
Most of the galaxies with counterparts have significant 
W1 and W2 flux measurements (98.83%). For W3 and 
W4, many of the flux measurements amount to upper 
limits (29.10% and 70.89%, respectively). The redshift 
distributions of galaxies with > 2<r W3 and W4 detec¬ 
tions are shown in Figure 1. We will provide flags to 
indicate which galaxies have detections, and which have 
upper limits. The usefulness of the upper limits will de¬ 
pend on the goals of the user. 

Galaxies are typically smaller than the WISE PSF, 
but the amount by which W?mpro underestimate the to¬ 
tal flux due to the spatial extent of the sources is not 
necessarily negligible. The Wl, W2, and W3 bands all 
have very similar PSFs with FWHM~ 6 — 7 arcsec, while 
the W4 PSF is larger (~ 12) arcsec. We investigate the 
missing light fraction for Wl, W2, and W3 by generat¬ 
ing with galfit (Peng et al. 2010) a series of simulated 
light distributions with 2-dimensional Sersic profiles con¬ 
volved with the Wl PSF models provided by Aniano 
et al. (2011), and inserted into empty sections of real 
Wl images 5 . Then, we use galfit to fit the PSF model 
to the simulated images to measure the PSF profile flux. 
We find that the difference between the PSF profile flux 
and the true, total flux is mostly a function of the effec¬ 
tive radius of the Sersic profile, and hardly depends on 
input flux, axis ratio, or Sersic index. For a typical size 
of 5 arcsec, we find that Ato varies from 15.63 mag for 
Sersic index n = 1 to 15.65 mag for Sersic index n = 4. 
Based on these simulations we find that the PSF profile 
magnitude underestimates the total flux by 

Am = 0.10 + 0.46 log(R e )+ 0.47 log(i? e ) 2 + 0.08 log(i? e ) 3 , 

(!) 

where R e is the effective radius in arcsec. This correction 
is used whenever the size is larger than R e = 0.5 arcsec. 
We propagate the uncertainty in the adopted radius into 
Delta m, which is in turn propagated into the flux un¬ 
certainty in flux. 

We use the r-band effective radii and apply a system- 

4 http://wise2.ipac.caltech.edu/docs/release/allwise/ 

5 http://www.astro.princeton.edu/~ganiano/Kernels/Ker_2012. 
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atic correction of a factor 1.5 ± 0.2 downward to to con¬ 
vert from optical size to near-IR size. This conversion is 
based on the analysis of (Vulcani et al. 2014), who con¬ 
sistently measured sizes for a large set of galaxies across 
the wavelength range 0.4 to 2.2 /im. These corrected sizes 
and Equation 1 are then used to correct the Wl, W2, and 
W3 fluxes. We use r-band effective radii are measured 
by Simard et al. (2011) when available, and otherwise 
deVRad or expRad from the SDSS catalog as appropri¬ 
ate. We prefer the Simard et al. (2011) measurements 
as their fitting methodology (two-dimensional light pro¬ 
files) is more similar to the methodology used by Vulcani 
et al. (2014) (the SDSS pipeline fits one-dimensional light 
profiles). We note that the difference between the two 
versions of our Wl/2/3 magnitude corrections is small: 
for objects for which both measurements are available 
the median difference is 0.03 mag, and the random scat¬ 
ter 0.14 mag. These uncertainties are small compared 
to the uncertainties in stellar masses and SFRs derived 
below. These final flux measurements, along with the 
applied corrections, are listed in the public catalog (see 
Section 4). The median correction is 0.25 mag, with a 
scatter of -0.07 (16%-ile) and +0.18 (84%-ile). Correc¬ 
tions for W4, for which the PSF is twice as wide, are not 
made. 

To account for possible systematic uncertainties of the 
WISE and SDSS photometric systems and further sys¬ 
tematic differences between the total flux measurements 
we adopt and propagate 0.1 mag uncertainties for all 
WISE fluxes. 

3. SED MODELING 
3.1. Method 

We use MAGPHYS to fit the photometric SED (da 
Cunha et al. 2008, 2012) 6 . The public version of 
MAGPHYS contains 50,000 stellar population template 
spectra (the optical photometric library) and 50,000 
PAH+dust emission template spectra (the infrared pho¬ 
tometric library). The stellar emission templates use the 
Bruzual & Chariot (2003) stellar population synthesis 
models and are generated for a wide range of star forma¬ 
tion histories parameterized as exponentially-declining 
models with superimposed random bursts of star forma¬ 
tion. 

The SEDs are computed by adding the individual spec¬ 
tra of all simple stellar populations, weighed by mass. 
The original MAGPHYS library was constructed for 
modeling IR-luminous, star-bursting galaxies. The orig¬ 
inal MAGPHYS library from the public version with 
50,000 stellar population template spectra was con¬ 
structed for modelling star-forming galaxies and we ex¬ 
tended the optical library to include 25,000 additional 
templates for more passive stellar populations. 

The infrared templates describe emission by dust. The 
total dust luminosity over 3 to 1000/xm has components 
of emission from PAHs and dust with a range of temper¬ 
atures. The model contains the ambient (diffuse) inter¬ 
stellar medium and star-forming regions (birth clouds). 
Because stars are born in dense molecular clouds which 
dissipate typically after 10 7 years, the SED of young stel¬ 
lar populations are attenuated by both dust in the birth 

6 http://www.iap.fr/inagphys 


clouds and dust in the ambient ISM, while the SED of 
older populations is only attenuated by dust in the dif¬ 
fuse ISM. The absorbed light is assumed to be re-emitted 
in the mid- and far-infrared, requiring conservation of en¬ 
ergy- 

MAGPHYS fits SEDs in the observed frame. There¬ 
fore, we generate libraries with model fluxes in the ob¬ 
served frame for a series of narrow redshift bins (Sz = 
0.0001) that span the range of our sample. The fit¬ 
ting process for each individual galaxy is then expedi- 
ated in two ways. First, we only consider optical tem¬ 
plates with similar ( g — i) m odei colors as the observed 
(g — i)data color. To be precise, we select templates with 
\(g-i)modei - ( g-i)data\ < 0.05 + cr[(g - i)data ], where a 
refers to the uncertainty. This method works for any two 
filters, and here we choose g — i. The selection is inclusive 
enough to avoid changes in the results: none of the mod¬ 
els eliminated from consideration have significant likeli¬ 
hood values. We check that we would get similar results 
if we fitted the whole model library. Second, we draw 
a random set of 1,000 infrared templates from the total 
set of 50,000 templates. Since our data do not sample 
the thermal peak, we cannot stringently constrain the 
dust temperature distribution, rendering full exploration 
of parameter space useless. As a test we compare the re¬ 
sults based on the full and reduced infrared libraries and 
find no systematic differences (smaller than 0.01 dex for 
both stellar mass and SFR) in the fitting results and no 
increases or decreases in the formal fitting uncertainties. 

In a small number of cases we find that MAGPHYS 
allows the presence of an unrealistically large amount of 
cold dust. In order to avoid this, we put a mild con¬ 
straint on the model 250/im luminosity based on the 
maximum observed 250/im-to-22/rm flux ratios (Chary 
& Elbaz 2001): Li' 250 fim[LQ/Hz\ < 85.7 x(Lv 22 /rm + 
cr(Li/ 22 ^m)) [ Lq/Hz\ , which is ~6 times the typical ra¬ 
tio. We will investigate the precision of our predicted 
total IR luminosities below, in § 3.3. 

Given our 10 flux constraints (SDSS ugriz; Wl—4; 
-bi+so/mi limit) x 2 is calculated for each model in the 
library. For each model parameter the posterior proba¬ 
bility distribution can be generated by taking the prior 
probability distribution (from the full library) and assign 
weights exp(—x 2 /2) to each of the models. In doing so 
we marginalize over all other model parameters and as¬ 
sume that the likelihood of a given model given the data 
is proportional to exp(— X 2 /2). We adopt the 50 per¬ 
centile value - that is, the median - as the best estimate. 
In practice, in the case of upper limits, we assumed zero 
flux with the upper limit as the error bar. 

Rest-frame fluxes are calculated using the observed 
fluxes and colors. Analogous to Holden et al. (2012) we 
derive for each redshift 2 bin a linear fit between, on the 
one hand, observed-frame magnitudes (m 0 b s 1 , rn 0 b s 2 ) of 
the templates in the filters straddling the desired rest- 
frame band, and, on the other hand, the rest-frame mag¬ 
nitude (mo) of the same templates: Too = m 0 bsi + A(z) 
x (rriobsi ~ Tn 0 bs 2 ) + B(z). Then, for each galaxy in 
the sample, we use the observed magnitudes and the val¬ 
ues of A(z) and B(z) to compute its rest-frame magni¬ 
tudes. We also adopt this technique to calculate rest- 
frame 12/mr and 22/mi luminosities. The basic result of 
the SED MAGPHYS modelling for each object, is a joint 
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PDF for the 16 parameters, which we then characterize 
by the median and the percentiles of the marginalized 
distributions. 

Fig. 2. — Rest-frame u — r vs. r — z distribution for our 
SDSS+WISE sample, color-coded according to the fraction of ob¬ 
jects that lie on the star-forming sequence or significantly below 
according to their full SED fit (see Section 4). The clearly distinct 
color-color distributions of star-forming and quiescent galaxies val¬ 
idate the assumption that a color-color diagram can be effectively 
used to distinguish the two types of galaxies, as is often done at 
high redshift. The labels correspond to the objects for which the 
SEDs are shown in Figures 3, 4, 5, and 6. 
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3.2. Sample SEDs 

The galaxies in our sample span a large range in stel¬ 
lar mass, star formation rates and colors (Figure 2). 
In Figures 3, 4, and 5 we show examples of SEDs and 
fitting results across the entire parameter space popu¬ 
lated by our sample. Figure 3 shows six massive galaxies 
(M* ~ 10 11 M 0 with little or no ongoing star formation 
(< 1 Mg yr _1 ). The characteristic red optical SEDs 
are accompanied by non-neglible mid-IR emission aris¬ 
ing from a combination of dust heated by evolved stars 
and traces of star formation. Figure 4 shows reddened 
galaxies with significant star formation activity, where 
the IR luminosity often exceeds the optical/NIR lumi¬ 
nosity. Figure 5 shows blue star-forming galaxies, with 
little or no extinction. 

To demonstrate the quality of the fits we show in Fig¬ 
ure 7 the chi-square distribution across parameter space 
(M*, SFR, colors). % 2 is the the goodness-of fit parame¬ 
ter for the best-fitting model, but not, formally speaking, 
the reduced chi-squared of the best-fitting model as dis¬ 
cussed by Smith et al. (2012). The generally low values 
of x 2 imply that our library of models is sufficiently ex¬ 
tensive. 


3.3. Verification of total IR Luminosities with Herschel 

MAGPHYS balances the energy absorbed in the 
UV/optical and the energy released in the IR, such that 
the inferred total IR luminosity can be expected to be 
correct if the optical/NIR SED is accurate. Still, since 
the WISE 12/un and 22/xm bands do not sample the ther¬ 
mal peak of the dust emission, it is important to test the 
robustness of the fitting results, in particular the star- 
formation rates and dust luminosities. 

Herschel-ATLAS provides PACS and SPIRE photom¬ 
etry over the wavelength range 100-500/nn. These obser¬ 
vations sample the total dust emission more directly than 
WISE. For the subsample of 285 galaxies that fall within 
the 16 square degree footprint of the publicly available 
part of Herschel-ATLAS (Eales et al. 2010; Ibar et al. 
2010; Pascale et al. 2011; Rigby et al. 2011; Smith et al. 
2011) we verify the precision and accuracy of our inferred 
total IR luminosity. 

Our SDSS+WISE-based models succesfully predict 
of total IR luminosities and star-formation rates from 
SDSS+WISE+Herschel photometry, as illustrated in 
Figure 6). Even though the marginalized probability 
distributions for the dust luminosity and star-formation 
rate are tightened when Herschel photometry is added, 
we find no systematic offset (0.00 dex) between the star- 
formation rates. Furthermore, the scatter (0.19 dex) is 
consistent with the formal confidence intervals (see Fig¬ 
ure 8). The uncertainties of SFR ~ 10~ 4 M Q /yr galaxies 
are large and dust luminosities are poorly constrained. 

3.4. The Public Catalog 

We provide two public catalogs 7 . They both contain 
all 858,365 galaxies from the SDSS spectroscopic galaxy 
sample with good redshift measurements as described in 
Section 2. Table 1 contains the input data (ID, redshift, 
fluxes, Galactic extinction). Table 2 contains the mod¬ 
eling results from MAGPHYS, V max (see below), rest- 
frame luminosities, and a set of flags. We recommend 
to use the MAGPHYS modeling results for objects with 
FLAG=1 (633,205 out of 858,365). These are all 2 < 0.2 
galaxies with reliable aperture corrections based on size 
measurements from Simard et al. (2011) ( FLAG-R = 1), 
good WISE photometry (FLAG-WI — 1 or 2), and 
good-quality SED fits (FLAG-CHI 2 = 1). This pri¬ 
mary sample contains 91.40% of all SDSS galaxies at 
2 < 0 . 2 . 

Our V max calculation only includes the maximum red¬ 
shift at which galaxies would be retained in the sam¬ 
ple. No low-redshift limit, related to the SDSS bright- 
magnitude limit, is taken into account. For our purposes 
of illustrating the distribution of SFR and stellar mass 
this does not matter given the large volume, but for in¬ 
dividual, low-redshift, high-luminosity galaxies the user 
should keep in mind that our V max cannot be directly 
used. 


4. DISCUSSION 

4.1. Comparison with Brinchmann et al. (2004) 

Stellar masses and SFRs from the MPA-JHU catalogs 
(Brinchmann et al. 2004) are widely used. The differ¬ 
ences between their and our stellar mass estimates are 

7 http://irfu.cea.fr/pisp/yu-yen.chang/sw.html 
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Fig. 3. — Example SEDs and fitting results of six massive, quiescent galaxies. The large sub-panels show the observed fluxes (red points), 
the best-fitting SEDs (black lines) and the corresponding unattenuated stellar SEDs (blue lines). The smaller panels show the marginalized 
posterior probability distributions of four models parameters: M*, SFR, dust luminosity (in solar units), and dust attenuation. The red 
lines show the median values, which we adopt as best-fitting values of the PDF. The short-hand IDs (Ml-6) correspond to the labels in 
Figure 2. 
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Fig. 4. — Example SEDs and fitting results of six star-forming galaxies , where the optical light is significantly dust-extincted. See Figure 
3 for an explanation of the panels, lines, and symbols. The short-hand IDs (Dl-6) correspond to the labels in Figure 2. 
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Fig. 5. — Example SEDs and fitting results of six blue, star-forming galaxies , with only modest dust extinction in the optical. See Figure 
3 for an explanation of the panels, lines, and symbols. The short-hand IDs (Sl-6) correspond to the labels in Figure 2. 
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Fig. 6.— Example SEDs and fitting results of six galaxies with Herschel-ATLAS photometry available (orange points) in addition to the 
SDSS &; WISE (red points). The blue and black lines represent fitting results to the SDSS + WISE data only (as in Figure 3), and the dark 
(light) gray areas represent the 68% (95%) confidence intervals for the far-IR dust SED predicted from these fits. The magenta and orange 
lines represent the best-fitting unattenuated and attenuated SEDs, respectively, when the Herschel-ATLAS SEDs are also simultaneously 
fitted. The smaller panels show the marginalized posterior probability distributions of four of the model parameters. In all cases, the 
SDSS+WISE fitting results predict the Herschel photometry well and all results inferred without Herschel data are consistent with the 
results inferred with Herschel data. However, the constraints on the dust luminosity and dust attenuation parameter are tightened once 
Herschel photometry is included. The short-hand IDs (Hl-6) correspond to the labels in Figure 2. 
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Fig. 7. — Quality of the SED fit, shown in dependence of SFR versus stellar mass (left; also see Figure 9 and Figure 11) and rest-frame 
u — r vs. r — z (right; also see Figure 2). The color coding is the median chi-square value. Boxcar smoothing in two dimensions considering 
the uncertainties of the parameters is applied to the color coding and opacity. For star-forming galaxies, the distribution is about uniform. 
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Fig. 8. — Top: SFR from the SED fits with only optical (SDSS) 
and near/mid-IR (WISE) data minus SFR from fits that also in¬ 
clude far-IR Herschel data vs. SFR from fits with Herschel data. 
The objects included here are those that overlap between our 
SDSS-selected sample and Herschel-ATLAS. Red error bars indi¬ 
cate 68% confidence intervals. Bottom: Same as above, but for dust 
luminosity. Our SDSS+WISE based fitting results are generally 
consistent with the SDSS+WISE+Herschel-based fitting results. 


the addition of WISE 3.4/im and 4.5/tm photometry as 
tracers of stellar mass, and the updated Galactic extinc¬ 
tion correction from Schlafly & Finkbeiner (2011). Fur¬ 
thermore, the simultaneous fitting of rnid-IR photometry 
provides a fundamentally different constraint on the dust 
properties of galaxies than optical photometry and emis¬ 
sion line ratios. In addition, dust attenuation laws are 
also different. In Brinchmann et al. (2004) the adopted 
slope for the attenuation curve is -0.7, following Chariot 
& Fall (2000). MAGPHYS has the same slope for diffuse 
dust, but a steeper slope (-1.3) for birth clouds. Despite 
these differences, offsets between the two sets of mass es¬ 
timates are small and insignificant (left-hand panel, Fig¬ 
ure 9). The reasons for this are that the same stellar 
population models and star-formation histories are used 
and that the impact of the WISE photometry is limited 
due to the relatively large uncertainties on the total flux 
measurements. 

Larger differences are seen for the SFR estimates. As 
opposed to the stellar mass estimates, the SFR estimates 
rely on wholly different tracers in the two cases. Selecting 
star-forming galaxies by their location in the color-color 
diagram (Figure 2), we find a median offset across the 
sample of 0.22 dex and a scatter of ~ 0.3 dex. This offset 
does not depend on mass or redshift, or even on whether 
the galaxies have significant 12/rm and 22yum detections, 
but strongly varies with SFR (Figure 9): for high-SFR 
galaxies our values are large compared to Brinchmann 
et al. (2004), for low-SFR galaxies our values are small. 
The anti-correlation between the changes in M* and SFR 
is the result of the underlying anti-correlation between 
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Fig. 9.— Comparison between the MPA-JHU values for M* and SFR and the values derived here (C15). Top left: median difference 
between the stellar mass estimates. These are always small (< 0.10 dex). Bottom left: scatter in the stellar mass estimates; generally 
< 0.15 dex. Top right: median difference in SFR. Here we only show star-forming galaxies selected the color-color classification given in 
the right-hand panel of Figure 2. Across the sample the offset is 0.22 dex. It is strongly SFR-dependent, and ranges from -0.3 dex to +0.3 
dex. Bottom right: scatter in the SFR estimates, typically 0.3 dex or a factor of two. The grey lines shows the star-forming sequence in 
Figure 11. 
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mass to light ratio and sSFR. 

4.2. 12ym and 22/im luminosities as SFR indicators 

Another consequence is that our 12/rm and 22^m SFR 
conversions are lower by, respectively, 0.22 dex and 0.10 
dex compared to the conversion based on a comparison 
with the Brinchmann et al. (2004) estimates, as carried 
out by Lee et al. (2013). 

In Figure 10 we show the correlations between 12^m 
and 22^m luminosities and our SFR estimates. The 
large downward scatter is due to mid-IR radiation that 
is not associated with star formation. In particular, the 
mid-IR luminosities of quiescent galaxies do not reflect 
star formation activity, but rather circumstellar dust and 
PAHs heated by evolved stars (e.g., Bressan et al. 2006). 
Even for many galaxies that are star-forming according 
to our definition - separated in color-color space (Fig¬ 
ure 2; SF galaxies: ( u — r) re st < 2.1 or ( u — r) res t < 
1.6 x (r — z) rest + 1.1) - a large fraction of the mid-IR 
luminosity is not attributed to star formation. Keeping 
this in mind, we derive the following conversions from 
mid-IR luminosity to SFR: 

log SFR/(M q yr- 1 ) = \ogL 12 /L Q - 9.18 (2) 

and 

log SFR/(Mq yr- 1 ) = log L 22 /L Q - 9.08, (3) 

These relationships are determined by calculating the 
median values of the SFR in bins of 0.25 dex wide bins 
in luminosity and performing a linear fit to these me¬ 
dian values, weighing by the inverse of the square root 
of the numer of objects in the bins. The upward scatter 
(84%-ile minus median) in SFR is 0.20 dex, while the 
downward scatter (median minus 16%-ile) is 0.60 dex, 
reflecting non-SF contributions. The scatter obviously 
depends on luminosity: above 10 10 L 0 the downward 
scatter is also 0.30 dex. 

The \2/m\ and 22ym values are 0.13 and 0.04 dex 
higher than the conversions provided by Jarrett et al. 
(2013), who based their SFRs on GALEX UV fluxes and 
the standard Spitzer/MIPS 24/zm calibration from Rieke 
et al. (2009). The relatively good agreement with the 
Jarrett et al. (2013) calibrations is encouraging, as their 
sample consists of 17 well-resolved, nearby galaxies for 
which measurement uncertainties are minimal. (In fact, 
the small difference is consistent with the random varia¬ 
tion expected in the average for small sample of 17 ob¬ 
jects.) Our calibration generalizes their result by extend¬ 
ing the dynamic range in stellar mass and star-formation 
rate by an order of magnitude upward and using a sample 
of several hundred thousand galaxies. 

Our 12/rm and 22/rm SFR conversions are lower by, re¬ 
spectively, 0.22 dex and 0.10 dex compared to Lee et al. 
(2013), who use the Brinchmann et al. (2004) SFRs. 8 
Part of this offset can be attributed to the systematic dif¬ 
ferences between our SFR estimates and those of Brinch¬ 
mann et al. (2004) (see Sec 4.1). Additional factors are 
the differences in sample selection and modeling tech¬ 
nique. We select star-forming galaxies on the basis of 

8 These differences include the appropriate correction from 
Salpeter to Chabrier IMF for the Lee et al. (2013) conversions. 


their optical colors, whereas Lee et al. (2013) select galax¬ 
ies with large 22/rm luminosities. Moreover, not all mid- 
IR radiation traces SF, which is taken into account in 
our modeling, whereas Lee et al. (2013) simply use total 
mid-IR luminosity as a SF tracer. 

4.3. A Fresh View on Bimodality 

Assessing bi-modality in the galaxy population re¬ 
quires a sample that is corrected for completeness and 
sufficient sensitivity in the star-formation tracer to sepa¬ 
rate star-forming and passive galaxies. We compare the 
magnitude in r band with the WISE bands at different 
redshift and find that the mass completeness limit of our 
sample is set by the SDSS r-band spectroscopic limit for 
all redshifts: the fraction of galaxies above the mass limit 
without WISE counterparts is less than 1%. The com¬ 
pleteness limit at a given redshift is found by identify¬ 
ing the most massive galaxies with magnitudes near the 
limit. As such we find that the mass completeness limit 
depends on redshift as log = 10.6+2.28 log(;z/0.1). 

The mass limit is such that actively SF-ing galaxies al¬ 
ways have significant 12^m detections; that is, if a galaxy 
is not detected at 12 /mi it must be a quiescent galaxy. 
In other words, the upper limit on SFR is useful is the 
sense that it allows us to distinguish between SF-ing and 
quiescent galaxies. Using the mass limit, a galaxy of a 
given mass is assigned a V ma x value assuming a survey 
area of 8,032 square degrees. These values are given in 
Table 2. 

In Figure 11 we show the SFR-M* distribution of the 
363,774 galaxies above our mass limit and with FLAG = 
1 (see Section 3.4). The sequence of star-forming galaxies 
(a.k.a. the Main Sequence) indicated in orange - is 
quantified by the median in 0.2-wide stellar mass bins. 
We fit these median values with a power law and find: 

log SFR/{M Q yr~ 1 ) = 0.80logM*/(10 10 M 0 ) - 0.23. 

( 4 ) 

We define the scatter in the Main Sequence as the 84%- 
50%-ile range, which is 0.39 dex. The downward scatter 
(50%-16%-ile range) is larger (0.64 dex), but this number 
is difficult to interpret due to the imperfect separation of 
star-forming and quiescent galaxies. This uncertainty 
is not (only) due to limited fidelity in our color-color 
classification, but also due to the natural variation in SF 
activity. 

For comparison, when we apply our fitting technique 
using the Brinchmann et al. (2004) values for stellar 
masses and SFRs, we find: 

log SFR/{M Q yr~ 1 ) = 0.75 logM*/(10 10 M 0 ) - 0.03. 

( 5 ) 

Our fit is slightly steeper, but more importantly, has a 
0.2 dex higher intercept at M* = 10 10 M 0 and a larger 
scatter (0.39 dex vs. 0.31 dex). The larger scatter can 
perhaps be attributed to the more direct measurement 
of highly obscured star formation when using mid-IR lu¬ 
minosities as a tracer as opposed to the optical emission 
lines. We note that the random uncertainties in our SFR 
estimates (typically, 0.18 dex) are much smaller than the 
scatter, indicating a large intrinsic scatter around the 
“main sequence” relation. This is in stark contrast with 
the measurements from Brinchmann et al. (2004): their 
typical uncertainty is 0.30 dex - the same as the scatter 
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Fig. 10. — SFR vs. rest-frame 12/zm and 22/xm luminosities for galaxies that have significant detections in both WISE bands. The data 
are binned and color-coded according to the fraction of star-forming galaxies as defined in the color-color diagram shown in Figure 2. For 
star-forming galaxies we produce a linear fit, as i ndicated by the solid lines. For comparison we show equivalent fits from Lee et al. (2013) 
and Jarrett et al. (2013). The SFRs are matched to Chabrier (2003) IMF. 
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Fig. 11.— Distribution of galaxies above our (redshift-dependent) 
mass completness limit in the SFR vs M* plane. The mass limit 
described in Section 4.3 is such that all galaxies on or near the 
the SF-ing sequence have significant detections at 12//rn. The star¬ 
forming sequence is shown in orange, where the dashed lines show 
the lfj scatter, as explained in Section 4.3. In order to visualize 
the SFR, bimodality we adopt the 1 -a upper limits on the SFR 
for galaxies where the SFR estimate is more uncertain than 0.39 
dex. Consequently, quiescent galaxies in this plot are placed at the 
upper limit of the true SFR. The gray scale represents the number 
of galaxies after weighing by 1 /Vmax ■ 


- implying zero intrinsic scatter. We conclude that the 


the star-forming sequence as revealed by the WISE pho¬ 
tometry is less tight than what is inferred from optical 
SDSS spectra. 

The bi-modality seen in Figure 11 in our inferred SFR- 
M* distribution can be directly compared with Figure 17 
and 24 from Brinchmann et al. (2004), where the SFRs 
are inferred from nebular emission lines strengths, and 
with Figure 15 from Salim et al. (2007), where the SFRs 
are inferred from UV luminosities. The color coding used 
in Figure 2 is based on the separation of star forming and 
non-star forming galaxies in Figure 11: non-star forming 
galaxies are those that lie below the sequence by 2 times 
the scatter or more. This confirms that the two-color 
diagram is a simple yet very effective tool to separate 
dusty from passive galaxies. 

Finally, it is of interest to explore the distribution of 
SF in relation to traditional tools to probe SF activity 
across the galaxy population, such as the color-mass di¬ 
agram, the BPT diagram (Baldwin et al. 1981), and the 
H5-D4000 relation (Kauffmann et al. 2003). In Figure 12 
we show the g — r vs. stellar mass distribution in three 
different ways: weighed by number density, weighed by 
SFR density, and weighed by stellar mass density. The 
first panel shows the well-known bimodality in the form 
of a red sequence and a blue cloud. We have also indi¬ 
cated the intermediate region often called the green val¬ 
ley, which has been argued to contain galaxies that are 
transitioning from blue to red (e.g., Schawinski et al. 
2007). However, as most recently shown by Taylor et al. 
(2015), many galaxies in this region are reddened be¬ 
cause of dust (also see, e.g., Conselice 2006) and not 
primarily due to a reduced level of SF. This is illustrated 
in the second panel, which shows that the galaxies in 
the green valley region - and not the more numerous, 
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Fig. 12. — Left: comoving number density distribution of galaxies in the stellar mass vs. rest-frame g — r plane. The grey scale reflects 
the numbers of galaxies in each bin after weighing with 1/Vmax • Middle: the grey scale reflects the SFR density. Right: the grey scale 
reflects the stellar mass density. Apparently, the green valley provides the most fertile ground for star formation as indicated by the green 
line. 


fainter blue galaxies - dominate the total SF budget of 
the present-day universe. On the other hand, these mas¬ 
sive, SF-ing galaxies do not dominate the stellar mass 
budget, as is shown in the third panel. The absence of 
a dominant population of massive, SF-ing galaxies and a 
general dearth of massive, disk-like galaxies (van der Wei 
et al. 2009) implies that their SFRs cannot be sustained 
for long (Schawinski et al. 2014). 

In Figure 13 we show the SFR distribution in the BPT 
and HS-D n ( 4000) diagrams in three different mass bins. 
In the left-hand panels, we show the BPT diagram at 
different stellar mass bins. For the low stellar mass sam¬ 
ple, most population are star-forming galaxies and most 
SF activities are also occur in SF region. For higher stel¬ 
lar mass bins, the population moves to composite galax¬ 
ies and even low-ionization nuclear emission-line region 
(LINER), but the SF activities still occur between star¬ 
forming and composite regions. 

In the right-hand panels we see, like Kauffmann et al. 
(2003), that the strength of the 4000Abreak, tracing 
evolved stellar populations (> 1 Gyr), is anti-correlated 
with HS, tracing the presence of younger stellar popula¬ 
tions (age < 1 Gyr). Not surprisingly, SF mostly occurs 
in galaxies with generally young populations. For low- 
mass galaxies the SF-ing population is representative of 
the general population, but for high-mass galaxies the 
SF-ing galaxies represent a tail of outliers, as most galax¬ 
ies are old and quiescent. Still, high-mass SF-ing galax¬ 
ies are older than their low-mass counterparts. This can 
be interpreted as evidence for an increasingly prominent 
bulge, while SF activity in the disk is similar to that in 
the disks of lower-mass galaxies as argued by Abramson 
et al. (2014). However, it could also signal a decline in 
SF activity as suggested by Schawinski et al. (2014). 

A detailed exploration of the processes that drive SF 
evolution is beyond the scope of this paper, but we hope 
that the community will take advantage of the mid-IR 
based SFR estimates published here to distinguish fur¬ 
ther between reddening by dust, age, and reduced SF 


activity. 

5. SUMMARY 

We revisit the measurement of stellar masses and star- 
formation rates for the SDSS spectroscopic galaxy sam¬ 
ple after adding 4-band photometry from WISE for the 
wavelength range 3 — 22^m. We benefit from this wave¬ 
length extension by adopting the latest, state-of-the-art 
SED modeling approach (MAGPHYS da Cunha et al. 
2008), which includes the self-consistent treatment of 
dust attenuation and emission and a wide range of star 
formation histories. The resulting SFR esimates are 
mostly based on PAH emission and thermal dust radi¬ 
ation, and, therefore, are complementary to the nebular 
emission line-derived estimates from (Brinchmann et al. 
2004). In Section 4.2 we provide new calibrations for the 
conversion of monochromatic 12/rni and 22/nn luminosi¬ 
ties into SFRs. 

The new M* and SFR estimates show the well- 
established bi-modality in SF activity. A first application 
is the verification of the cruder color-color separation of 
SF-ing and quiescent galaxies that is often used at higher 
redshifts. We confirm that the U — V-V — J (or, alterna¬ 
tively, u — r-r — z, as in Figure 11) color-color separation 
has a 89% success rate in identifying SF-ing and quies¬ 
cent galaxies. 

We also fit the relation between SFR and M* (the SF- 
ing sequence, or SF-ing Main Sequence). Our sequence 
is slightly steeper than previous studies based on SDSS 
galaxies, and has a higher normalization and scatter. We 
attribute these differences to our new SFR estimates that 
now include dust emission and also to some extent the 
uncertainty in establishing a purely star-forming sample 
(i.e. without contamination from quiescent galaxies). L* 
falaxies contribute most to the cosmic SF density, and 
these typically have optical colors that are in between 
those of the traditional red squence and blue cloud (Fig¬ 
ure 12) as a result of dust attenuation, not because of 
reduced star-formation activity. 
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Fig. 13.— Number density (contours) and SFR density distribution (grey scale) in the BPT diagram (left, only galaxies with detected 
emission lines) and Ui5-D n (4000) (right). From top to bottom we show 3 different mass bins. 


As a service to the community we make our 
SDSS+WISE matched photometry catalog as well as the 
SED fitting results publicly available. 
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TABLE 1 
Input Catalog 


Column Name Format Unit Column Description 

NYU-VAGC catalog index 
J2000 R.A. [deg] from NYU-VAGC (r-band) 

J2000 Dec. [deg] from NYU-VAGC (r-band) 

Redshift from the NYU-VAGC spectroscopic catalog 


ID 

LONG 

- 

RA 

DOUBLE 

deg 

DEC 

DOUBLE 

deg 

REDSHIFT 

DOUBLE 

- 

PLATE 

LONG 

- 

MJD 

LONG 

- 

FIBERID 

LONG 

- 

DESIGNATION 

STRING 

- 

FLUXO.U 

DOUBLE 

Jy 

FLUX0.U.E 

DOUBLE 

Jy 

FLUXO.G 

DOUBLE 

Jy 

FLUX0.G.E 

DOUBLE 

Jy 

FLUX0.R 

DOUBLE 

Jy 

FLUXCLR-E 

DOUBLE 

Jy 

FLUXO-I 

DOUBLE 

Jy 

FLUX0J.E 

DOUBLE 

Jy 

FLUX0.Z 

DOUBLE 

Jy 

FLUX0.Z.E 

DOUBLE 

Jy 

FLUX0.W1 

DOUBLE 

Jy 

FLUX0.W1.E 

DOUBLE 

Jy 

FLUX0.W2 

DOUBLE 

Jy 

FLUX0.W2.E 

DOUBLE 

Jy 

FLUX0.W3 

DOUBLE 

Jy 

FLUX0.W3.E 

DOUBLE 

Jy 

FLUX0.W4 

DOUBLE 

Jy 

FLUX0.W4.E 

DOUBLE 

Jy 

FLUX.U 

DOUBLE 

Jy 

FLUX_U_E 

DOUBLE 

Jy 

FLUX.G 

DOUBLE 

Jy 

FLUX_G_E 

DOUBLE 

Jy 

FLUX.R 

DOUBLE 

Jy 

FLUX.R.E 

DOUBLE 

Jy 

FLUX.I 

DOUBLE 

Jy 

FLUXJJ5 

DOUBLE 

Jy 

FLUX.Z 

DOUBLE 

Jy 

FLUX_Z_E 

DOUBLE 

Jy 

FLUX.W1 

DOUBLE 

Jy 

FLUX.W1.E 

DOUBLE 

Jy 

FLUX.W2 

DOUBLE 

Jy 

FLUX_W2_E 

DOUBLE 

Jy 

FLUX.W3 

DOUBLE 

Jy 

FLUX.W3.E 

DOUBLE 

Jy 

FLUX.W4 

DOUBLE 

Jy 

FLUX_W4_E 

DOUBLE 

Jy 

EXTIN.U 

DOUBLE 

mag 

EXTIN.G 

DOUBLE 

mag 

EXTIN.R 

DOUBLE 

mag 

EXTIN.I 

DOUBLE 

mag 

EXTIN.Z 

DOUBLE 

mag 


SDSS plate from NYU-VAGC 

SDSS mjd from NYU-VAGC 

SDSS fiberid from NYU-VAGC 

ALLWISE designation 

22 -band Flux before corrections 

Uncertainty of 22 -band flux before corrections 

g-band flux before corrections 

Uncertainty of g-band flux before corrections 

r-band flux before corrections 

Uncertainty of r-band flux before corrections 

2 -band flux before corrections 

Uncertainty of 2 -band flux before corrections 

z-band flux before corrections 

Uncertainty of z-band flux before corrections 

W 1-band flux before corrections 

Uncertainty Wl-band of flux before corrections 

W 2-band flux before corrections 

Uncertainty W2-band of flux before corrections 

W 3-band flux before corrections 

Uncertainty W3-band of flux before corrections 

W4-band flux before corrections 

Uncertainty W 4-band of flux before corrections 

u-band Flux after corrections 

Uncertainty of 22 -band flux after corrections 

g-band Flux after corrections 

Uncertainty of g-band flux after corrections 

r-band Flux after corrections 

Uncertainty of r-band flux after corrections 

z-band Flux after corrections 

Uncertainty of 2 -band flux after corrections 

z-band Flux after corrections 

Uncertainty of z-band flux after corrections 

W 1-band Flux after corrections 

Uncertainty W1-band of flux after corrections 

W 2-band Flux after corrections 

Uncertainty W 2-band of flux after corrections 

W 3-band Flux after corrections 

Uncertainty W 3-band of flux after corrections 

W 4-band Flux after corrections 

Uncertainty W4-band of flux after corrections 

Galactic extinction correction for 22 -band 

Galactic extinction correction for g-band 

Galactic extinction correction for r-band 

Galactic extinction correction for 2 -band 

Galactic extinction correction for z-band 
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TABLE 2 
Output Catalog 


Column Name 

Format 

Unit 

Column Description 

LMASS.2.5 

FLOAT 

log Mg 

log stellar mass (2.5th percentile) 

LMASS.16 

FLOAT 

log Mg 

log stellar mass (16th percentile) 

LMASS.50 

FLOAT 

log Mg 

log stellar mass (50th percentile) 

LMASS.84 

FLOAT 

log Mg 

log stellar mass (84th percentile) 

LMASS.97.5 

FLOAT 

log Mg 

log stellar mass (97.5th percentile) 

LSFR.2.5 

FLOAT 

log Mg/jr 

log SFR. (2.5th percentile) 

LSFR.16 

FLOAT 

log Mg/jr 

log SFR (16th percentile) 

LSFR.50 

FLOAT 

log Mg/jr 

log SFR (50th percentile) 

LSFR.84 

FLOAT 

l°gMg/V 

log SFR (84th percentile) 

LSFR.97.5 

FLOAT 

log Mg/jr 

log SFR. (97.5th percentile) 

LSSFR.2.5 

FLOAT 

logl /yr 

log specific SFR (2.5th percentile) 

LSSFR.16 

FLOAT 

log 1 /yr 

log specific SFR (16th percentile) 

LSSFR.50 

FLOAT 

log 1 /yr 

log specific SFR (50th percentile) 

LSSFR.84 

FLOAT 

log 1 /yr 

log specific SFR (84th percentile) 

LSSFR.97.5 

FLOAT 

logl /yr 

log specific SFR (97.5th percentile) 

LDUST.2.5 

FLOAT 

log Lq 

log dust luminosity (2.5th percentile) 

LDUST.16 

FLOAT 

log Lq 

log dust luminosity (16th percentile) 

LDUST.50 

FLOAT 

log Lq 

log dust luminosity (50th percentile) 

LDUST.84 

FLOAT 

log Lq 

log dust luminosity (84th percentile) 

LDUST.97.5 

FLOAT 

log Lq 

log dust luminosity (97.5th percentile 

MU.2.5 

FLOAT 

~ 

dust attenuation parameter in da Cunha et al. (2008) (2.5th percentile) 

MU_16 

FLOAT 

~ 

dust attenuation parameter in da Cunha et al. (2008) (16th percentile) 

MU-50 

FLOAT 

~ 

dust attenuation parameter in da Cunha et al. (2008) (50th percentile) 

MU.84 

FLOAT 

~ 

dust attenuation parameter in da Cunha et al. (2008) (84th percentile) 

MU-97.5 

FLOAT 

~ 

dust attenuation parameter in da Cunha et al. (2008) (97.5th percentile) 

TAUV-2.5 

FLOAT 

- 

dust attenuation parameter in da Cunha et al. (2008) (2.5th percentile) 

TAUV-16 

FLOAT 

- 

dust attenuation parameter in da Cunha et al. (2008) (16th percentile) 

TAUV.50 

FLOAT 

- 

dust attenuation parameter in da Cunha et al. (2008) (50th percentile) 

TAUV-84 

FLOAT 

- 

dust attenuation parameter in da Cunha et al. (2008) (84th percentile) 

TAUV-97.5 

FLOAT 

- 

dust attenuation parameter in da Cunha et al. (2008) (97.5th percentile) 

V.MAX 

DOUBLE 

Mpc 3 

Maximum volume for LMASS-50 

LREST-U 

DOUBLE 

log Lq 

Rest-frame zz-band Luminosity 

LREST.G 

DOUBLE 

log Lq 

Rest-frame g-band Luminosity 

LREST.R 

DOUBLE 

log Lq 

Rest-frame r-band Luminosity 

LREST.I 

DOUBLE 

log Lq 

Rest-frame z-band Luminosity 

LREST-Z 

DOUBLE 

log Lq 

Rest-frame z-band Luminosity 

LREST.W1 

DOUBLE 

log Lq 

Rest-frame IF1-band Luminosity 

LREST.W2 

DOUBLE 

log Lq 

Rest-frame W 2-band Luminosity 

LREST.W3 

DOUBLE 

log Lq 

Rest-frame W 3-band Luminosity 

LREST-W4 

DOUBLE 

log Lq 

Rest-frame W 4-band Luminosity 

FLAG-R 

INT 

- 

1: Simard radius; 2: deVaucouleurs radius; 3: exponential radius; 0: no radius 

FLAG-W 

INT 

- 

1: single optical matched in WISE; 2: > 1 counterparts within 6”; 0: not matched 

FLAG-W1 

INT 

- 

1: detected W1 (SNR>2); 2: upper limit on Wl; 0: no W1 data 

FLAG.W2 

INT 

- 

1: detected W2 (SNR>2); 2: upper limit on W2; 0: no W2 data 

FLAG.W3 

INT 

- 

1: detected W3 (SNR>2); 2: upper limit on W3; 0: no W3 data 

FLAG-W4 

INT 

- 

1: detected W4 (SNR>2); 2: upper limit on W4; 0: no W4 data 

FLAG-CHI2 

INT 

- 

1: x 2 < 3 for best-fit model; 0: y 2 > 3 for best-fit model 

FLAG 

INT 

- 

1: good fit (FLAG R=1; FLAG-W?=1 or 2; FLAG-CHI2=1; 2 < 0.2); 0: others 











